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Abstract 



A general unifying framework for integrable soliton-like systems on time scales is 
introduced. The i?-matrix formalism is applied to the algebra of 5-differential operators 



in terms of which one can construct infinite hierarchy of commuting vector fields. The 



theory is illustrated by two infinite-field integrable hierarchies on time scales which 
are difference counterparts of KP and mKP. The difference counterparts of AKNS and 
Kaup-Broer soliton systems are constructed as related finite-field restrictions. 

m 

1 Introduction 

qq ■ Integrable systems are widely investigated in (1 + 1) dimensions, where one of the dimensions 
stands for the time evolution variable and the other one stands for the space variable. The 
space variable is usually considered on continuous intervals, or both on integer values and 
on R [lj or on K q intervals [2j [3] . In order to embed the study of integrable systems into a 
more general unifying framework, one of the possible approaches is to construct the integrable 
systems on time scales. Here the space variable is considered on any time scale where R, KL, 
K g are special cases. The first step in this direction was taken in where the Gelfand- 
Dickey approach E] was extended in order to construct integrable nonlinear evolutionary 
equations on any time scale. Another unifying approach is to formulate different types of 
discrete dynamics on R. Some contribution in this direction was made recently in [7]. 

The main goal of this work is to present a theory for the systematic construction of (1 + 1)- 
dimensional integrable systems on time scales in the frame of the i?-matrix formalism. By 
an integrable system, we mean such a system which has an infinite-hierarchy of mutually 
commuting symmetries. The i?-matrix formalism is one of the most effective and systematic 
methods of constructing integrable systems [SI [9] . This formalism originated from the pio- 
neering article [3] by Gelfand and Dickey, who constructed the soliton systems of KdV type. 
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The crucial point of the .R-matrix formalism is that the construction of integrable systems 
proceeds from the Lax equations on appropriate Lie algebras [ED, [9] . The simplest .R-matrices 
can be constructed by a decomposition of a given Lie algebra into two Lie subalgebras. We 
refer to [HI El H] for abstract formalism of classical i?-matrices on Lie algebras. 

This paper is organized as follows: In the next section, we give a brief review of the time 
scale calculus. In the third section, we define the (^-differentiation operator and formulate the 
Leibniz rule for this operator. We introduce the Lie algebra as an algebra of 5-differential 
operators equipped with the commutator, decompose it into two Lie subalgebras and con- 
struct the simplest i?-matrix on this algebra. We present the appropriate Lax operators for 
infinite-field cases and the admissible finite-field restrictions generating consistent Lax hier- 
archies. In T = R case, or in the continuous limit of some special time scales, we observe 
that the algebra of 5-differential operators turns out to be the algebra of pseudo-differential 
operators. Next, we formulate and prove the property of the algebra of ^-differential op- 
erators. This property allows us to obtain natural constraints which are fulfilled by finite 
field restrictions. Therefore, the source of the constraints, obtained in the Burgers equations 
and KdV hierarchy on time scales in [3j, is established. We end up this section with the 
construction of the recursion operators by means of the method presented in [10J. In the 
fourth section, we illustrate two infinite-field integrable hierarchies on time scales which are 
difference counterparts of Kadomtsev-Petviashvili (KP) and modified Kadomtsev-Petviashvili 
(mKP) hierarchies. In the last section, we present finite-field restrictions which are difference 
counterparts of Ablowitz-Kaup-Newell-Segur (AKNS) and Kaup-Broer (KB) hierarchies with 
their recursion operators. 



2 Preliminaries 



In this section, we give a brief introduction to the concept of time scale. We refer to [TTl [T2"] 
for the basic definitions and general theory of time scale. What we mean by a time scale T, is 
an arbitrary nonempty closed subset of real numbers. The time scale calculus was introduced 
by Aulbach and Hilger [T3l [H] in order to unify all possible intervals on the real line K, like 
continuous (whole) M, discrete Z, and g-discrete K q (K q = q z U {0} = {q k : k G Z} U {0}, 
where q ^ 1 is a fixed real number) intervals. For the definition of the derivative in time 
scales, we use forward and backward jump operators which are defined as follows. 

Definition 2.1 For ieT, the forward jump operator a : T — > T is defined by 

a(x) = inf {y G T : y > x}, (2.1) 

while the backward jump operator p : T — > T is defined by 

p(x) = sup {y G T : y < x}. (2.2) 

We set in addition cr(maxT) = maxT if there exists a finite maxT, and p(minT) = minT if 
there exists a finite minT. 

The jump operators a and p allow the classification of points in a time scale in the following 
way: x is called right dense, right scattered, left dense, left scattered, dense and isolated if 



a(x) = x, a(x) > x, p(x) = x, p{x) < x, a(x) = p(x) = x and p(x) < x < cr(x), respectively. 
Moreover, we define the graininess functions p, v : T — > [0, oo) as follows 

p(x) = a(x) — x, v(x) = x — p(x), for all x G T. (2.3) 

In literature, T K denotes a set consisting of T except for a possible left-scattered maximal 
point while T K stands for a set of points of T except for a possible right-scattered minimal 
point. 

Definition 2.2 Let f : T — > R be a function on a time scale T. For x G T K , delta derivative 
of f , denoted by Af, is defined as 

A/(«)=]Jm MzZM , seT , (2.4) 

u>/ii/e /or x G T K , V -derivative of f , denoted by V f , is defined as 

V/ W = Um /(s) - / . ( f )) , S6 T, (2.5) 

provided that the limits exist. A function f : T — ► R is said to &e A-smooth (V -smooth) if it 
is infinitely A- differentiable (V -differentiable). 

Remark 2.3 Let / : T — > R 6e A- differentiable on T K . If x is right- scattered, then the 
definition ( 12. 4ft torns out to 6e 

A / (X) = /KX) . ) ~ /(X) , 
u>/ii/e i/x is right-dense, ( 12.41) implies that 

A f(x) = lim s G T. 

Similarly, let f : T — > R be V-differentiable on T K . If x is left- scattered, then the definition 
(12.51) turns out to be 

(V(x) 

u>/ii/e i/x is left-dense, (12.51) yields as 

rrt, \ i- /(g) ~ /( g ) c tt 

Vj(x) = hm , s G T. 

s— >a; X — S 

In order to be more precise, we present A and V derivatives for some special time scales. If 
T = R, then A- and V-derivatives become ordinary derivatives, i.e. 

A/On) = V/(x) = ^M. 



If T = HZ, then 
If T = Kg, then 

A J x = —7 7\ and V/ x = — 3- , 

(q — l)x (1 — q L )x 

for all x 7^ 0, and 

A/(0) = V/(0) = fim /(S) ~ /(Q) , s G Kg, 
provided that this limit exists. 

As an important property of A-differentiation on T, we give the product rule. If /, g : T — > R 
are A-differentiable functions at x G T K , then their product is also A-differentiable and the 
following Lebniz-like rule hold 

A(fg)(x) = g(x)Af(x) + f(a(x))Ag(x) 

= f(x)Ag(x)+g(a(x))Af(x). 1 ' ' 

Besides, if / is A-smooth function, then 

f(a(x))=f(x)+v(x)Af(x). (2.7) 

If x G T is right-dense, then fi(x) = and the relation (12.71) is trivial. 

Definition 2.4 A time scale T is regular if both of the following two conditions are satisfied: 

(i) a(p(x)) = x for all x G T, 

(ii) p(a(x)) = x for all x G T. 

Set = minT if there exists a finite minT, and set 2* = —00 otherwise. Also set x* = maxT 
if there exists a finite maxT, and set x* = 00 otherwise. 

Proposition 2.5 jl] A time scale is regular if and only if the following two conditions hold: 

(i) the point = minT is right dense and the point x* = maxT is left-dense; 

(ii) each point ofT \ is either two-sided dense or two-sided scattered. 



In particular R, fiZ (H 7^ 0) and K q are regular time scales, as are [0, 1] and [—1,0] U {1/k 
k G N} U {k/(k + 1) : k G N} U [1, 2]. 



Throughout this work, let T be a regular time scale. By A, we denote the delta-differentiation 
operator which assigns each A-differentiable function / : T — ► R to its delta-derivative A(/), 
defined by 

[A(f))(x) = Af(x), for xeT K . (2.8) 
The shift operator E is defined by the formula 

(Ef){x) = f(a(x)), xeT. (2.9) 

The inverse E~ x is defined by 

(E- l f)(x) = f(a-\x)) = f(p(x)), (2.10) 

for all x ET. Note that E~ x exists only in the case of regular time scales and that in general 
E and E~ x do not commute with A and V operators. 

Proposition 2.6 [15j Let T be a regular time scale. 

(i) If f : T — > R is a A-smooth function on T K , then f is V -smooth and for all x G T K; 

V/(x) = E~ l Af{x). (2.11) 

(it) If f : T — > R is a V -smooth function on T K , then f is A-smooth and for all x G T K , 

Af(x) = EVf(x). (2.12) 

Thus the properties of A- and V-smoothness for functions on regular time scales are equiva- 
lent. 

In some special cases, by properly introducing the deformation parameter, it is possible to 
consider a continuous limit of a time scale. For instance, the continuous limit of fflL is the 
whole real line R, i.e. 

T = KL -5=^ T = R; (2-13) 
and the continuous limit of K q is the closed half line R + U 0, thus 

T = K q T = R + U0. (2.14) 

For more about the calculus on time scales we refer the readers to [HI [T2"] . 



3 Algebra of ^-differential operators 
3.1 Basic notions 

In this section, we deal with the algebra of ^-differential operators defined on a regular time 
scale T. We denote the delta differentiation operator by 5 instead of A, for convenience in the 



operational relations. The operator Sf which is a composition of 5 and /, where / : T — > R, 
is introduced as follows 

6f:=Af + E(f)6, V/. (3.1) 
Note that, the definition (13.1 j) is consistent with the Lebniz-like rule on time scales (12.61) . 

Theorem 3.1 The Leibniz rule on time scales for the operator 5 is given as follows, 
(i) For n ^ : 

n 

5 n f = J2 Yl (A ik + 1 EA ik E...A h EA h )fd k , (3.2) 

k=0 ii+i2+...+ik+i=n—k 

where i 1 ^ for all 7 = 1,2, .., k + 1. Here the formula includes all possible strings 
containing n — k times A and k times E. 

(ii) For n < : 

00 

S n f=Y^ (-l) k+n (E~ ik + n+1 AE- lk +"A...E- i2 AE- ll )f5- k , (3.3) 

k=-n u+«2+---+«fc+n+i=fc 

where i 1 > for all 7 = 1, 2, .., /c + n + 1 > 0. i/ere £/ie formula includes all possible 
strings containing k + n + 1 tzmes and k + n times A. 

The above theorem is a straightforward consequence of definition (13.11) . Note that S~ 1 f has 
the form of the formal series 

00 

r 1 / = ^2(-l) k ((E- 1 A) h B- 1 )f6- h -\ (3.4) 

fc=0 

which was previously given in @], in terms of V. Thus (13.31) is the appropriate generalization 

of (J33D. 

3.2 Classical i?-matrix formalism 

In order to construct integrable hierarchies of mutually commuting vector fields on time 
scales, we deal with a systematic method, so-called the classical R-matrix formalism [HIEIQ], 
presented in the following scheme. 

Let Q be an algebra, with some associative multiplication operation, over a commutative field 
K of complex or real numbers, based on an additional bilinear product given by a Lie bracket 
[',']'■ G G, which is skew-symmetric and satisfies the Jacobi identity. 



Definition 3.2 A linear map R : Q — > Q such that the bracket 

[a, o]r := [Ra, b] + [a, Rb], 
is a second Lie bracket on G , is called the classical R-matrix. 



(3.5) 



Skew-symmetry of (13.51) is obvious. When one checks the Jacobi identity of (13. 5p . it can be 
clearly deduced that a sufficient condition for R to be a classical i?-matrix is 

[Ra, Kb] - R[a, b] R + a[a, b] = 0, (3.6) 

where a G K, called the Yang-Baxter equation YB(a). There are only two relevant cases of 
YB(ct), namely a^O and a = 0, as Yang-Baxter equations for a ^ are equivalent and can 
be reparametrized. 

Additionally, assume that the Lie bracket is a derivation of multiplication in Q, i.e. the 
relation 

[a, be] = b[a, c] + [a, b]c a,b,c G Q (3.7) 

holds. If the Lie bracket is given by the commutator, i.e. [a, b] = ab — be, the condition (13.71) 
is satisfied automatically, since Q is associative. 

Proposition 3.3 Let Q be a Lie algebra fulfilling all the above assumptions and R be the 
classical R-matrix satisfying the Yang-Baxter equation, YB(a) . Then the power functions L n 
on Q , L G Q and n G Z +; generate the so-called Lax hierarchy 

^- = [R(L n ),L], (3.8) 

of pairwise commuting vector fields on Q. Here, t n 's are related evolution parameters. We 
additionally assume that R commutes with derivatives with respect to these evolution param- 
eters. 



Proof. It is clear that the power functions on Q are well defined. Then 

[L tm )t n ~ (L tn ) tm = [RL m ,L] tn - [RL n ,L] tm 

= [(RL m )t n - (RL n )t m ,L] + [RL m , [RL n , L]] - [RL n , [RL m , L]] 
= i(RL m ) tn - (RL n ) tm + [RL m , RL n ],L]. 

Hence, the vector fields (I3.8P mutually commute if the so-called zero- curvature (or Zakharov- 
Shabat) equations 

{RL m ) tn - {RL n ) tm + [RL m , RL n ] = 0, 

are satisfied. From ( 13.81) and by the Leibniz rule ( 13.71) we have that (L m ) tn = [RL n ,L m ]. 
Using Yang-Baxter equation for R and the fact that R commutes with dt n , we deduce 

R(L m )t n -R(L n ) tm + [RL m ,RL n ] = 

= R[RL n , L m ] - R[RL m , L n ] + [RL m , RL n ] 

= [RL m , RL n ] - R[L m , L n ] R = -a[L m , L n ] = 0. 

Hence, the vector fields pairwise commute. □ 

In practice the powers of Lax operators in (13. 8p are fractional. Notice that, the Yang-Baxter 
equation is a sufficient condition for mutual commutation of vector fields (13.81) . but not nec- 
essary. Thus choosing an algebra Q properly, the Lax hierarchy yields abstract integrable 
systems. In practice, the element L of Q must be appropriately chosen, in such a way that 
the evolution systems (13.81) are consistent on the subspace of Q. 



3.3 Classical i?-matrix on time-scales 



We introduce the algebra Q as an algebra of formal Laurent series of (pseudo-) (^-differential 
operators equipped with the commutator, and define its decomposition such as: 

o = g >k ® Q<k = {J2 ® iJ2 u i( x ) 5i }> ( 3 - 9 ) 

C^-k i<k 

where Ui : T — > K are A-smooth functions. The subspaces G^k, Q<k are closed Lie subalgebras 
of Q only if k — 0, 1. Thus, we define the classical i?-matrix in the following form 

R:=l(P >k -P <k ) k = 0,1, (3.10) 

where P^fc and P<fc are the projections onto G^k and respectively. Since the classical 
P-matrices (I3.10P are defined through the projections onto Lie subalgebras, they satisfy the 
Yang-Baxter equation (13.61) for a — j. 

Let L G Q be given in the form 

L = u N 5 N + un-^' 1 + . . . + mS + u + M_i<5 -1 + . . . , (3.11) 

where U{ are dynamical fields depending additionally on the evolution parameters t n . Thus, 
the Lax hierarchy ( 13. 8ft , based on (13.101) and in general generated by fractional powers of L, 
turns out to be 

§ L n = [( Lf )^> L ]=-[( L ^< k > L ] k = ^ 1 nGZ +- ( 3 ' 12 ) 

Proposition 13.31 implies that the hierarchy (13.121) is infinite hierarchy of mutually commuting 
vector fields and represents (1 + l)-dimensional integrable differential- difference systems on a 
time scale T, including the time variables t n and space variable iGl 

Analyzing (13.121) for L given by (13. lip , in the case of k — 0, one finds that (w^v)* = and 
(■UAT_i)i = /i(. . .) (see also Remark 14.11) . Similarly for k = 1, we have (mat)* = /x(. . .) (see 
also Remark 14.21) . Hence, the appropriate Lax operators, yielding consistent Lax hierarchies 
(13.121) . are in the following form: 

k = 0: L = c N 5 N + u N . 1 S N ~ 1 + .., + utS 1 +u + u^5- 1 + ... (3.13) 
k = l: L = u N 5 N + un^- 1 + . . . + ux8 l + u + v-rf- 1 + . . . , (3.14) 

where cn is a time-independent field and fields un-i,un are time-independent for dense x G T, 
as at these points ji = 0. This is the reason why they are distinguished by a tylde mark. 

Nevertheless, we are interested in finite-field integrable systems on time-scales. Thus, in order 
to work with a finite number of fields, we should impose some restrictions on (I3.13P and (13.141) 
in such a way that the commutator on the right-hand side of the Lax equation ( I3.12p does not 
produce terms not contained in the left-hand side of the Lax equation. To be more precise, 
the left- and right-hand of ( 13.121) span the same subspace of Q. From this purpose, in the 
case of k = 0, one finds the general admissible form of finite-field Lax operator given by 

L = c N S N + mjv-i^" 1 + . . . + u x 8 + M + ^s S ~ X( P» ( 3 - 15 ) 



with further restriction 

L = c N S N + mtv-i^ 1 + • • • + ui5 + uq. (3.16) 
In the case of k = 1, the general admissible Lax operator has the form 

L = u N 5 N + un-iS*' 1 + . . . + u t S + u + 5~ x u- x + V, ( 3 - 17 ) 

s 

and further restrictions are 

L = u N S N + mtv-i^" 1 + • • • + u x 8 + M + <J -1 u_i (3.18) 
L = u N 5 N + un^S*- 1 + . . . + uiS + u (3.19) 
L = u N 5 N + UN-tS"- 1 + . . . + ui5. (3.20) 

In the above Lax operators cjy is a time-independent field for all x G T and un-i^un are 
time-independent at dense points from a time scale. We assume also that the sum J2 S ^ s 
finite. 

In general, for an arbitrary regular time scale T, the Lax hierarchies (13.121) represent hier- 
archies of soliton-like integrable difference systems. For instance, when T = YiL or K g , the 
hierarchies f)3.12p are those of lattice and g-deformed (-like) (discrete) soliton systems, respec- 
tively. In particular, for the case of T = M, i.e. the continuous time scale on the whole M, 
the Lax hierarchies are those of field soliton systems. In some cases, field soliton systems can 
also be obtained from the continuous limit of integrable systems on time scales (see (12.131) 
and (J2HD). 

In the continuous time scale, the algebra of 5-differential operators (13.91) turns out to be the 
algebra of pseudo-differential operators 

{Y,u l (x)d l }®{Y,u l (x)d 1 }, (3.21) 

ijjfc i<k 

where d is such that du = d x u + ud = u x + ud. The above decomposition is valid only 
if k = 0, 1 and 2. Thus, in the general theory of integrable systems on time scales, we 
loose one case in contrast to the ordinary soliton systems constructed by means of pseudo- 
differential operators. This follows from the fact that, for k = 2, (13.91) does not decompose 
into Lie subalgebras for an arbitrary time scale. For appropriate Lax operators, finite field 
restrictions and more information about the algebra of pseudo-differential operators, we refer 
the reader to [161 El El H]. Note that the fields i/j s and (p s in (13.151) and (13.171) are special 
dynamical fields in the case of the algebra of pseudo-differential operators. They are the so- 
called source terms, as tp s and (p s are eigenfunctions and adjoint-eigenfunctions, respectively, 
of the Lax hierarchy ( 13.121) [T?] . 

It turns out that there are constraints between dynamical fields of the admissible finite-field 
Lax restrictions ( 13. 15ft - (13. 201) fulfilling (13. 12ft . We give these constraints in the following 
theorem, which is a consequence of the property of the algebra of 5-differential operators. 
This property is illustrated in the following lemma. 
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Lemma 3.4 Consider the equality 

r 

5 r F = J2Ci5 r -\ r>0. (3.22) 

i=0 

Then the following relation 

r 

^(-/i) 1 ^ = F (3.23) 

i=0 

is valid. 

Proof. We make use of induction. Assume that (13.231) holds for r. Then 

r r r+1 

S r+1 F = S r( EF ^ § + S r AF = ^ ^-i+l + ^ ^.jr-i = ^ Q^+l"*. (3.24) 

i=0 i=0 i=0 

By the assumption we have X^I=o( — /-O'A = -^-^ an d Si=o( — A*) 1 -^* = AF. Hence 

r+1 r r 

^(-/i) J Q = ^(-/i)^ 1 ^, + ^2{-fj)% = —fiAF + EF = F. (3.25) 

i=0 i=0 i=0 

□ 

Let us explain the source of Lemma I3.4L Consider the equality 

A = ^(^ = 0, (3.26) 

where the sum is finite, and A is purely 5-differential operator. We expand A with respect to 
the shift operator £: £u = E(u)S. From the relation (12.71) we have 

£ = 1 + ^5. (3.27) 

The equality from Lemma 13.41 is trivially satisfied for dense i£T, since in this case /i = 0. 
Thus, it is enough to consider remaining points in a time scale so assume that [i ^ 0. Hence, 
from (13.271) . we have the formula 

5 = iM- 1 S-fM- 1 . (3.28) 
Thus, using (13. 28f) the relation (I3.26f) can be rewritten as 

A = J2a' i £ i = 0. (3.29) 

i 

Obviously, it must hold for terms of all orders. The equality for the zero-order terms, i.e. 
a' = 0, can be simply obtained by replacing 5 with — /i -1 in (13.261) . The same substitution in 
(13.221) allows us to find 

r 

(-^ F = J2C t (-^r r+ \ (3.30) 

i=0 



which is equivalent to (13.231) . 

The above procedure can be extended also to operators A that are not purely 5-differential 
and contain finitely many terms with 5~ 2 , .... As an illustration consider the equality 

r-l 

[A5 r , ip5~ V] = Yl °i sr ~ 1 ' i + CrS- 1 ^ + i)5- l C r . (3.31) 

The above equality is well-formulated since it follows immediately from the definition and the 
property of the 8 operator. Replacing 5 with — the commutator vanishes, and we have 

r-l 

= ^a(-/u)- r ' +1+i + C r {-fi)(p + fjL)C r (3.32) 

i=0 

r-l 

Y^i-vYCi + (-/u) r (d.^ + V>C r ) = 0. (3.33) 
Straightforward consequence of such a behavior of 5-differential operators is the next theorem. 



Theorem 3.5 



(i) The case k = 0. The constraint between dynamical fields of (13. 15ft . generating Lax 
hierarchy (I3.12p . has the form 

^ ^N^ duN-i + ^ydUj _ ^ ^ d ^ sip ^ = 

dt n dt n dt n 

N _ 2 S (3-34) 

(-//) jv_1 m JV -i + ^(-/u) l Mi - fj,^2ip a (p a = a n , 

i=0 s 

where n 6 Z + and a n zs a time-independent function. 

(ii) The case k = 1. T/ie constraint between dynamical fields of ( 13.171) . generating ( 13.121) . 
has the form 

(3-35) 



i=-l 



where n E Z + and a n is a time-independent function. 



Proof. We already know that Lax operators (13.151) and (13.171) generate consistent Lax hi- 
erarchies ( 13.121) . Thus, the right-hand side of ( 13.121) can be represented in the form of L tn . 
Replacing 5 with — \T X in ( 13.121) . we have 

LtJ s =-u-i= [(L n hk,L}L 1= 0. (3.36) 



Hence, the constraints (13.341) and (13.351) follow. 



□ 



The above theorem can be generalized to further restrictions. As a consequence, the con- 
straints ( I3.34p or ( 13.35!) with fixed common value of all (In-, £11*6 valid for the whole Lax 
hierarchy (I3.12p . 



3.4 Recursion operators 

One of the characteristic features of integrable systems possessing infinite-hierarchy of mu- 
tually commuting symmetries is the existence of a recursion operator [THl H|- A recursion 
operator of a given system, is an operator of such property that when it acts on one sym- 
metry of the system considered, it produces another symmetry. Giirses et al. [TDJ presented 
a general and very efficient method of constructing recursion operators for Lax hierarchies. 
Among others, the authors illustrated the method by applying it to finite-field reductions of 
the KP hierarchy. In [19] the method was applied to the reductions of modified KP hierarchy 
as well as to the lattice systems. Our further considerations are based on the scheme from 
[ID] and PS]. 

The recursion operator $ has the following property: 

HLt n ) = L tn+N , neZ + , 

and hence it allows reconstruction of the whole hierarchy (13.121) when applied to the first 
(N — 1) symmetries. 

Lemma 3.6 

(i) The case k = 0. Let the Lax operator be given in the general form (13. 15)) . Then, the 
recursion operator of the related Lax hierarchy can be constructed solving 

L tn+N =L tn L+[R,L) (3.37) 

with the remainder in the form 

R = on^- 1 + ■ ■ ■ + a Q + a-i, s rV s , (3.38) 

where N is the highest order of L. 

(ii) The case k — 1. Similarly for the Lax operator ( 13. 171) the recursion operator can be 
constructed from (13.371) with 

R = a N 5 N + ■ ■ ■ + a + a_i jS rV s - (3.39) 



-L- '-J 



Proof. Consider the case k = 0. Then for f 1 3 . 1 5 j) we have 
(L^)^o = {{L%) >0 L) >0 + ((L^) <0 L)^ 

s 

where [£\ a5 l ]o = and R is given by (I3.38p . Similarly for k = 1, we have 
(L^)^ = ((Lt)^L)^ + ((^)<i^)^i 

s 

= (L*) >X L + R, 

where R has the form (13. 39ft . Thus, in both cases (13.371) follows from (I3.12p . Hence we can 
extract the recursion operator from (13. 37ft . □ 

Note that in general, recursion operators on time scales are non-local. This means that they 
contain non-local terms with A -1 being formal inverse of A operator. However, such recursion 
operators acting on an appropriate domain produce only local hierarchies. 

4 Infinite-field integrable systems on time scales 

4.1 Difference KP, k = 0: 

Consider the following infinite field Lax operator 

L = 6 + u + ^2ui6-\ (4.1) 

which generates the Lax hierarchy ( 13. 121) as the difference counterpart of the Kadomtsev- 
Petviashvili (KP) hierarchy. 

For (L)^ = 5 + u , the first flow is given by 

^ = J2(-l) h+1 Ui- k J2 (E-^AE-^A...E^AE^)u ( 4 - 2 ) 

1 fc=0 ji+j2 + -+jk+i=i 

+ fiAu i+ i + Aui + UiU Vz > 0, 
where j 7 > for all 7 ^ 1. 

For (L 2 )>o = 8 2 + £5 + T], where 

£:=Eu + u r] := Au + u 2 + u x + Eui, (4.3) 



one calculates the second flow 

fiA(E + l)u 2 + fj,A(Aui + uiUq + u%E~ u Q ) 



du 
dt 2 

dui 
dt 2 



i-l 

J2 (E-^AE-^A . . . E~ j2 AE- jl )£ 
k=-i h+j2+-+j k +2=i+i (4.4) 

i-l 

+ ^(-l^+V-fc (E~ jk+1 AE~ jk A . . . E~ j2 AE~ jl )t] 

k=0 jl+j2+—+jk+l=i 

+ A 2 Ui + (EA + AE)u i+l + fj,A(E + l)u i+2 + £(A Ui + Eu i+l ) + 7]u h 
where j 7 > for all 7 ^ 1. 

The simplest case in (2 + 1) dimensions: We rewrite the first two members of the first flow 
by setting u = w and t\ = y and the first member of the second flow by setting t 2 = t. Since 
E and A do not commute, note that in the calculations up to the last step, we use E — 1 
instead of /iA, to avoid confusion. 

w y = {E-l)u x , (4.5) 
u 1>y = (E-l)u 2 + A Ul + Ul (l- E- l )(w), (4.6) 
w t = {E 2 -l)u 2 + {E -\){Au 1 + u l w + u 1 E- 1 (w)) (4.7) 

Applying E + 1 to (14. 6 p from the left yields: 

(E 2 - l)u 2 = (E + l)u 1>y — (E + 1)A«! -{E- l)u x {\ - E- l )w. (4.8) 

Applying (E — 1) to (14.7ft from the left and substituting ( 14.51) and (14.81) into the new derived 
equation we finally obtain the (2 + l)-dimensional one- field system of the form 

fiAw t = (E+ i)w yy - 2 Am y + 2fiA(ww y ). (4.9) 

which does not have a continuous counterpart. For the case of T = hZ, one can show that 
( 14.91) is equivalent to the (2 + l)-dimensional Toda lattice system. 

The difference analogue of one-field continuous KP equation is too complicated to write it 
down explicitly. 

Remark 4.1 Here we want to illustrate the behavior of Uq in all symmetries of the difference 
KP hierarchy. Let (L n ) <0 = v^d^ 1 , then by the right-hand of the Lax equation (I3.12p . we 

obtain the first members of all flows 

f^A<,<*>. (4.10) 



Thus Uq is time-independent for dense x G T since /i = 0. Hence when T = R, uq appears to 
be a constant. 
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In T = M. case, or in the continuous limit of some special time scales, with Uq = 0, the Lax 
operator ( 14.1 ft turns out to be a Laurent series of pseudo-differential operators 

L = d + s ^u l d- i . (4.11) 
Moreover, the first flow (14. 2p turns out to be exactly the first flow of the KP system 

= V * > 1 ( 4 - 12 ) 



while the second flow ( 14.41) becomes exactly the second flow of the KP system 



d i_1 / — 1\ 

'"' 'u i ) 2x + 2{u i+1 ) x + 2j2(-±) k+1 ^ k Jui-k{ui) k x Vi^l. (4.13) 



dt 2 

fc=i 



4.2 Difference mKP, = 1: 

Consider the Lax operator of the form 

L = u-i5 + ^Mi5 _< (4.14) 

which generates the difference counterpart of the modified Kadomstsev-Petviashvili (mKP) 
hierarchy. 

Then, (L)^i = U-\8 implies the first flow 



dti 



^= ^(-i) fc +V- fe (E- jk +*AE- j ^A...E- j *AE-^)u„ 1 ( 415 ) 

1 fc=-i h+h-\ — hifc +2 =«+i 

+ u_iEu i+ i + u^iAui Wi ^ 0, 

where j 7 > 0, 7 = 1, 2, . . . , k + 2. 

Next, for (£ 2 )^i = £<5 2 + i]5, where 

£ := U-iEu-i, 7] := u^iAu-i + u^iEuq + uqu^i, (4.16) 

we have the second flow as follows 
du-i 



dt 2 
dui 
dt 2 



£(EAu + E l {ui)) + {iu^Auq - «i£ £ - uliAuo 

i-l 

{-l) k+3 Ui- k J2 (E~ jk+3 AE~ jk+2 A . . . AE- jl )£ 
k=-2 j'i+i2+...+ifc+3=«+2 (4.17) 

+ (-l) k+2 Ui- k (E~ jk+2 AE~ jk+1 A . . . AE- jl )r] 

k=-x ii+i2+---+jfe+2=*+i 
+ i 2 {A 2 Ui + (EA + AE)u i+1 + E 2 u i+2 ) + r](A Ul + Eu i+1 ), 

where i ^ and j 7 > for all 7 ^ 1. 



Remark 4.2 Similarly in order to illustrate the behavior ofu_i in all symmetries of the dif- 
ference mKP hierarchy let us consider (£ n )<i = v^8~ % . Then we obtain the first members 

of all flows 

^ = MU_xAt;W (4.18) 

Thus u-x is time-independent for dense x G T. Hence when T = R ; u_i appears to be a 
constant. 

In T = K. case, or in the continuous limit of some special time scales, with u-\ = 1, the Lax 
operator ( 14. 14ft turns out to be the pseudo-differential operator 

L = d + ^2uid-\ (4.19) 

Furthermore, the system of equations ( 14.151) is exactly the first flow of the mKP system 

= u i!X , \/i > 0, (4.20) 
while the second flow ( 14.17P turns out to be the second flow of the mKP system 



-jr = (ui) 2x + 2(u i+ i) x + 2u (ui) x + 2u u i+1 

0X2 

+ 2j2(-l) k+1 ( l \ l+1 ^ k (u ) kx Vi ^ 0. 

1 — n \ / 



(4.21) 



5 Finite-field integrable systems on time scales 
5.1 Difference AKNS, k = 0: 

Let the Lax operator (13.151) for N = 1 and c\ — 1 is of the form 

L = 5 + u + ip5- 1 <f. (5.1) 
The constraint (13.341) between fields, with a n = 0, becomes 

u = flip (p. (5.2) 

For (L)j> — 8 + u, one finds the first flow 

^ = /.A(^-V), 

^ = ^ + A^, (5.3) 
^f- = -wp + A£TV 



Eliminating field u by (15.21) we have 



dip , 2 . , 
— = fxip ip + Aip, 
at i 

= + A£~V 

at i 



(5.4) 



Next we calculate (L 2 )^ = 5 2 + £5 + i] where 

£:=(E+l)u, r]:= Au + u 2 + ^E(ip)+ipE- l (^). (5.5) 
Thus, the second flow takes the form 

(III 

— = fiA [A(ipE~ l ((p)) + iPE- l (u(f) + uipE^ip] - fiA(E + l^E' 1 AE' 1 ^) 
ati 

^r- = tpn + £Aip + A 2 ip (5.6) 
dt 2 

^ = - i pr ] + AE-\^)-(AE- 1 ) 2 i f. 

By the use of the constraint ( 15. 2p . the second flow can be written as 

^ = ip(AfMip(p + (^<p) 2 + (fE(ip) + ipE~ l ( v )) + (E + l)iiipipAip + A 2 ip, 
dt 2 

%f = -<p(AfiiJ>tp + (iiipip) 2 + (pE(if>) + *pE~ l ((p)) + AE- l (ip(E + l)fiil><p) - {AE- l fip. 
ai2 

(5.7) 

In order to obtain the recursion operator one finds that for the Lax operator (15.11) the appro- 
priate reminder ( 13.381) has the form 

R = A- 1 (/i- 1 ^) - ^rV- (5-8) 

Then, (13.371) implies the following recursion formula as 

u-fi- 1 <pE ipE- 1 \ /u\ 

ip + ipA- 1 ^- 1 A + u + ipA^ip ipA^ip ^ ( 5 - 9 ) 

-ip A- 1 fi' 1 -cpEA^cp u — AE^ 1 — (pEA~ l ip J \p) t 

valid for isolated points x G T, i.e. when fi ^ 0. For dense points one must use its reduction 
by constraint (15.21) : 







V 


1 -( 




tn+1 



ip\ (A + u + fitpip + 2ip A~ V /j,ip 2 + 2ipA~ 1 ijj \ ftp 

<J ~{ -if(E + l)A- l ip u — AE' 1 — (p(E + l)A~ l ip) \<p J . 

where u is given by (15. 2p . 



(5.10) 



In T = R case, or in the continuous limit of some special time scales, with the choice u = 0, 
the Lax operator (15.11) takes the form L = d + ipd~ 1 ip. Then, the continuous limits of (15.31) 
and (15.6j) respectively, imply that the first flow is the translational symmetry 



— = ip 
dti 

dip 

dti x 

and the first non-trivial equation from the hierarchy is the AKNS equation 

dip o 

-j— = i^xx + 2ip <f, 
dt 2 

d<£ „ 2, 

TT = "Vase - 2y? 1p- 

dt 2 

For that special case the recursion formula (I5.10I) is of the following form: 

ip\ fd x + 2ipd- 1 <p 2ipd~ 1 'ip \ ftp 



ip) u+i \ -2ipd x V ~d x - 2<pd x V/ \V 



(5.11) 



(5.12) 



(5.13) 



5.2 Difference Kaup-Broer, k = 1: 

^From the admissible finite field restrictions ( 13. 17ft . we consider the following simplest Lax 
operator 

L = u5 + v + 5~ 1 w. (5.14) 
The constraint (I3.35p . with a n = 1, implies 

u = 1 + fxv — fj?w. (5.15) 



Then, for = u5, the first flow is given as 

du 

— = fiuAv, 
dti 



dv 

— =uAv + ^AE-\uw), (5.16) 
dti 

dw . „ i / „ v 

— = AE (uw). 
dti 



By the constraint ( 15.151) one can rewrite the first flow as 

dv 

—— = (fiv — fi 2 w)Av + fiAE (w(/iv — /i 2 w)), 
dti 

^ = AE' 1 (fivw - fi 2 w 2 ) . 



(5.17) 



Next, we calculate (£ 2 )^i = £<5 2 + r]S, where 

£ := uEu, r\ := uAu + uEv + vu, 



(5.18) 



that yields the second flow 
du 

—— = nuA(E~ l + l)uw + fiuAv 2 + fiuA(uAv), 
dt 2 

dv 

— = £(A 2 v + Aw) + fiAE-\w7]) + E^AE-^wg) + rjAv, (5.19) 
at 2 

^ = -AE^AE-Hwt) + AE-\wr 1 ). 
dt 2 

One can rewrite the above system reducing it by the constraint, but the final equation has 
complicated form. 

For the Lax operator (15.141) the appropriate reminder f!3.39j) is given by 

R = uA~ l {iiu)- l u t J - v tn - A- l w tn . (5.20) 

Hence, from (13.371) we have the following, valid when fj, ^ 0, recursion formula 

Ruu uE fj,u \ I u\ 

v + uA (l + E-^u lv\ , (5.21) 

w —AE~ l u-\-v — iiwJ \w J 

where 

Run = E(v) - yT x u + pLuA{v)A~ x {pM)~ x 

R vil = A(v) +w + MA(t»)A" 1 (/iM)" 1 + (1 - E~ l )uwA- l (nu)- 1 (5.22) 
R W u = AE^uwA^tfiuy 1 . 

Its reduction by the constraint ( 15.151) is 







• 


1 ■( 


w 


tn+1 





v\ _ (v + uA + R vii fx (I + E - Rvufj 2 \ (v 
wj. \ w + R W u^ -AE~ 1 u + v - aw - R wii fi 2 J \w 



(5.23) 



with u given by (15.151) . 



In the case of T = R, or in the continuous limit of some special time scales, with the choice 
u = 1, the Lax operator ( 15.141) takes the form L = d + v + d~ 1 w. Then the similar continuous 
analogue allows us to obtain the first flow 

dv 
dti 
dw 

drr Wx > 

and the first non-trivial equation from the hierarchy is the Kaup-Broer equation 



(5.24) 



dv 

-jr = v 2x + 2w x + 2vv x , 
dt 2 

dw 

-jr = ~w 2x + 2(vw) x . 
at 2 

For such special cases, the recursion formula (15.231) turns out to be 

v\ = fd x + v + v x d~ 1 2 \ / v 
w J \ w + dxwd' 1 —d x + v J \w 



(5.25) 



(5.26) 
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